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ABSTRACT 

Using the spectroscopic catalogue of the Sloan Digital Sky Survey Data Release 10, 
we have explored the abundance of satellites around a sample of 254 massive (10 11 < 
M* <2x10 u Mq) local (z < 0.025) galaxies. We have divided our sample into four mor¬ 
phological groups (E, SO, Sa, Sb/c). We find that the number of satellites with M* >10 9 Mq 
and R<300 kpc depends drastically on the morphology of the central galaxy. The average 
number of satellites per galaxy host (Afsat/ATiost) down to a mass ratio of 1:100 is: 4.5 + 0.3 
for E hosts, 2.6 ± 0.2 for SO, 1.5 + 0.1 for Sa and 1.2 + 0.2 for Sb/c. The amount of stellar 
mass enclosed by the satellites around massive E-type galaxies is a factor of 2, 4 and 5 larger 
than the mass in the satellites of SO, Sa and Sb/c-types, respectively. If these satellites would 
eventually infall into the host galaxies, for all the morphological types, the merger channel 
will be largely dominated by satellites with a mass ratio satellite-host /r >0.1. The fact that 
massive elliptical galaxies have a significant larger number of satellites than massive spirals 
could point out that elliptical galaxies inhabit heavier dark matter haloes than equally massive 
galaxies with later morphological types. If this hypothesis is correct, the dark matter haloes 
of late-type spiral galaxies are a factor of ~2-3 more efficient on producing galaxies with the 
same stellar mass than those dark matter haloes of early-type galaxies. 

Key words: galaxies: abundances - galaxies: elliptical and lenticular, cD - galaxies: evolu¬ 
tion - galaxies: formation - galaxies: luminosity function, mass function - galaxies: spiral 


1 INTRODUCTION 

Galaxy mergers have been raised in the last years as the most likely 
channel of size and mass growth of massive (M* >10'‘M©) galax¬ 
ies through cosmic time. Numerous observational and theoretical 
studies support this mode of growth, a mechanism that has in¬ 
creased the size and mass of the massive galaxies during the last 
~ 11 Gyr. In this scenario, the ancestors of the present-day most 
massive galaxies created the bulk of t heir mass in a short but very 


inten s e starburst event at z>2 (e.g . Keres et al 


20091: lOser et al Il2010l: iRicciardelli et al]l20l(J:Twuvts et al.ll20ld: 


Dekel et al] 


Boumaud et al] 1201 if) having, in that first evolutionary stage, a 


structure more compact. Later, a progressive process of merg¬ 
ers with satellites produ ced the envelopes that we see today sur¬ 
rounding these galaxi es dKhochfar & Silkl I 2 OO 6 I : lOser et alJl20ld : 
IFeldmann et all |201 II) . Many works support the above scheme, 
finding evidences for a continuous size evolution of the massive 
galaxies since z~3 fe.g. iTruiillo et alJ[2007l : iBuitrago et alJl2008h 
mainly produced by the formation of the outer most regions (e.g. 
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afl hoiol : 

In addition, other au¬ 
thors have found that the average velocity dispersion of the massive 
galaxies have decreased mildly since z~ 2 as expected theoreticall y 
from the galaxy merger scenario (e.g ICenarro & Truiillol l2009h . 
Other observations point out that this size evolution of the mas- 
sive galaxies does n ot depend on the age of the stellar population 
(iTruiillo et alfeoi lb nor on their intrinsic sizes dDiaz-Garcfa et al] 
12013 ). This suggests a growth engine external to the galaxy prop¬ 
erties. The absence of a significant number of relic galaxies in the 
nearby Universe also favours a merging scenario (e.g. ITruiillo et al.l 
I2009LI20 1 4 iTavlor et al .11201(1) . 

The above merging scenario can be alternatively probed mea¬ 
suring the satellite abunda nces around massive galaxies as cosmic 
time flows (e.g. lNewman et afl2012l) . There are many works that 
have studied in detail the properties of the satellite galaxies over 
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2014]). In particular. 


Marmol-Queralto et al. ( 20121) found that the fraction of massive 


galaxies with satellites of a given mass ratio (1:100 up to ; = 
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1 and 1:10 up to z = 2) have remained constant with time. 
A behaviour which is in qualitative good agreement with semi- 
analytical predicti ons based on the A cold dark matter (ACDM) 
model ( lOuilis & Truiillol2012l ). However, the semi-analytical mod¬ 
els over-predict the fraction of massive galaxies with satellites 
down to 1:100 mass ratio by a factor of ~2. 

Parameters such as the abundance of satellites, their distribu¬ 
tion or their intrinsic properties are intimately bound up with their 
host merger histories. These properties are thus, closely related to 
the underlying cosmology and they can be used to establish use¬ 
ful constrains to the models. The colours and structural proper¬ 
ties of the host galaxies can be modified by gravitational interac¬ 
tions with their satellites. The main goal of this work is to analyse 
the relation between the abundance of satellites and the host mor¬ 
phology in a sample of nearby massive (~10 u Mq) galaxies. We 
segregate our sample of galaxies in four groups which are repre¬ 
sentative of different structural configurations. Our morphological 
classification identifies visually the groups E, SO, Sa and Sb/c to 
explore the correlation between the number of satellites and the 
morphology of the h ost. Our approach differs from previous stud- 
ies based on colours dChenll2008l : lGuo et al ]|2011bl: Wane & White! 
l2012h . or those based on more general grouping (e.g. e arly-, late- or 
spheroid- disc like) of the massive galaxies (e.g. Guo et al]|201 lbl : 


g. inio 

120121 ) . 


Marmol-Oueralto et al]|2012l:lNierenberg et ali 

Our samples of massive galaxies could be used in future works 
to test the ACDM predictions about the number of satellites sur- 
roundin g the most massive galaxies in the present-day Universe 
(see e.g l lchenll2008l:lLiu et alJl201ll : IWang & Whitell2012l) accord¬ 
ing to their morphology. Our study also allows us to explore which 
is the most likely merging channel of present-day massive galax¬ 
ies, i.e. which type of satellites contribute most to the mass increase 
of their host galaxies in case they eventually merge with its main 
object. Consequently, this local stu dy, along with other works at 
higher z (see e.g. iFerreras et all2014l) . allows us to explore whether 
the merging channel has changed with time and its dependence with 
the host morphology. 

This paper is structured as follows. In Section [2] we describe 
our sample of hosts and satellite galaxies, their completeness and 
their stellar mass estimates. Section [3] explains the satellite selec¬ 
tion criteria and the methods used to clean our sample from back¬ 
ground and clustering contamination. Our results concerning the 
satellite abundances of the distinct samples are presented in Sec¬ 
tion [4] Section [5] discusses the main results of this paper and fi¬ 
nally our work is summarized in Section [6] Hereafter, we assume 
a cosmology with O m = 0.3, Qa = 0.7 and Ho = 70 km s _1 Mpc~'. 


dure to estimate the stellar masses of the galaxies of the catalogue. 
Then, we describe the selection of our host galaxies, the catalogue 
of potential satellites and finally, we study the completeness of our 
satellite population. 


2.1 The stellar mass estimation 

One of the goals of this work is to analyse the abundance of satel¬ 
lites as a function of the mass_ratio satellite-host. To estimate the 
stellar masses , we use the [Bell et alj l i2003l) ’s recipe, assuming a 
Kr oupal (i200l i) initial mass function (IMF). We take the ‘model- 
Mag’ g- and r-band magnitudes from the ‘specphoto’ SDSS DR10 
catalogue once they have been corrected from Galactic extinction 
dSchlegel et al.lll998i) . Then, we estimate the mass-to-light (MIL) 
ratio from the rest-frame g—r colour, being the M/L ratio estimated 
in the r band as follows: 

log (M/L), = a r + b r (g - r) - 0.15, (1) 

where a, = -0.306 and b r = 1.097 are the coefficients applied 
for determining the M/L ratio in the r band and 0.15 is subtracted 
to match the results to a Kroupa IMF. 

Using these computed (M/L),, we can directly estimate the 
stellar masses using the next relationship: 

log(M/M 0 ) = log (M/L) r - 0.4(M r - M 0 , r ), (2) 

where M r is the absolute magnitude of the galaxy and 
M 0 , r =4.68 the absolute magnitude of the Sun in the SDSS r band. 
Given that our study is focused on objects at very low redshift, we 
do not apply A'-corrections to the above g and r magnitudes since 
it will not affect significantly our mass estimates. In fact, the ex¬ 
pected values for ATcorrections at z < 0.025 are typically below ~1 
per cent relative to calibrati on errors found for g and r filters in the 
photometry of SDSS DR10 JPadmanabhan et alj||2008l ). 

To test how reliable our stellar mass estimates based on 
colours are, we h a ve co mpared our stellar masses wi t h thos e from 
iNair & Abrahaml J2010h based on iKauffmann et alj <1 20031) . This 
comparison can be done only for a subset of the galaxies in our 
cat alogue. As the res ult of this comparison we find a bias of 0.2 
dex INair & Abrahaml (being the [2010, smaller) with a typical un¬ 
certainty of 0.1 dex among the two stellar mass estimators. On what 
follows, we take that uncertainty as representative of our error es¬ 
timation of the stellar mass. The above bias is not surprising tak¬ 
ing into account the different methodologies and stellar population 
models used in both estimates of the stellar mass. 


2 THE DATA 

In this paper, we use the ‘specphoto’ spectroscopic catalogue 
of Sloan Digital Sky Survey (SDSS) Data Release 10 (DR10; 
Ahn et al.l 120141) to explore the abundance of satellites around a 
sample of massive galaxies in the nearby Universe (z < 0.025). The 
spectroscopic completeness of this catalogue is 90 per cent down 
to r = 17.7 mag. The catalogue includes a total of 1507954 Bary- 
onic Oscillation Spectroscopic Survey (BOSS) spectra comprising 
927844 galaxy spectra, 182009 quasar spectra and 159327 stellar 
spectra selected over 6373.2 deg 2 . We select those objects labelled 
as ‘GALAXY’ within the data set ‘specphoto’. This subset only has 
galaxies where the ‘SpecObj’ is a ‘sciencePrimary’ object, and the 
BEST PhotoObj is a PRIMARY. 

We structure this section as follows. First, we show the proce- 


2.2 The sample of host galaxies 

Using the available data, we estimate the stellar mass of our 
galaxies as we explained above (Section l2.lt . To build our sam¬ 
ple of massive host galaxies, we select only the galaxies with 
M* >10 u Mq. We limit our sample to galaxies with z < 0.025 (i.e. 
at a distance <100 Mpc). The average apparent r band magnitude 
of our host galaxies is r ~13 mag. This implies, taking into account 
the spectroscopic completeness limit of our catalogue, that we can 
identify potential satellites with stellar masses 100 times less mas¬ 
sive than their hosts. We have discarded massive hosts whose M/L 
(as computed by their g — r colour) were unreasonable (i.e. out of 
the interval 0.1 < (M/L), < 4) and those ones with large pho¬ 
tometric errors. The percentage of massive galaxies discarded by 
these reasons were 5 per cent. We also check visually our sample 
of massive galaxies to reject objects wrongly labelled as galaxies 
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Figure 1. Distribution of the normalized stellar mass for the host of each 
morphological type. The upper left-hand panel corresponds to the massive 
E type, the upper right to SO type, the lower left to Sa and the lower right to 
the Sb/c massive galaxies. The non-stripped area illustrates the distribution 
of massive galaxies chosen (10 11 <M*/Mq <2x10 i 1 ) for this work. The 
Njiost label in the panels indicates the number of host found within that 
stellar mass interval. 



Figure 2. Four examples of massive galaxies in our sample and their sur¬ 
roundings. The dashed circles enclose a region of 300 kpc in radius. The 
small circles indicate the position of the satellites. The upper right-hand 
panels are a zoom-in of the host galaxies to illustrate their morphologies. 


as well as those ones in clear interaction with another galaxy. Also, 
to avoid incompleteness in the number of satellites due to area cov¬ 
erage, 27 massive galaxies were additionally discarded due to their 
proximity to the edge of the SDSS DR10 spectroscopic footprint. 

Finally, the mass distribution of the different morphological 
types was checked. Our morphological segregation is explained in 
the next subsection. The distribution of the mass of the hosts is 
illustrated in the Fig. [T] and shows a high-mass-tail for the massive 
ellipticals which is not so prominent for the rest of galaxy types. To 
avoid any potential bias caused by this different mass distribution of 
the hosts, we established an upper mass limit of 2xlO n Mg. Thus, 
our final host sample is composed of 254 massive galaxies. 

The main goal of this work is to study the local abundance 
of satellites as a function of the morphological type of the massive 
hosts. Unlike other studies whose separation is based on colours, 
we do our own visual classification based on the Flubble classifica¬ 
tion using the SDSS DR10 images (see next subsection). Within the 
mass range (10 11 <M*/Mq <2x 10 u ) we find 83 E-type, 38 S0- 
type, 81 Sa-type and 52 Sb/c-type galaxies. Fig.[2]illustrates some 
of our host g alaxies. Our visua l classi fication is compared with the 
one done bv lNair & Abraham] d2010l) for 89 massive galaxies we 
have in common. We find an agreement of 84 per cent. 


2.2.1 Sorting the galaxy hosts into morphological classes 

In this work we have segregated our galaxies into four different 
morphological types (E, SO, Sa and Sb/c). The physical motiva¬ 
tion for sorting into these four categories is related to the expected 
strong connection between the evolutionary path followed by the 
galaxies and their detailed morphologies. This connection leaves 
their imprints on the relation between the specific angular momen¬ 
tum of the g alaxies, at a fixed stellar ma ss, and the galaxy morphol¬ 
ogy (see e.g. lRomanowskv & Falll2()T3) . This link can also be seen 
in the different shape of the outermost regions of the galaxies de¬ 


pending on the global m orphological type (e.g. Poh len & Truiillol 
l2006l : lErwin et al.l2008l) . If as expected, the merging activity of the 
galaxies is connected to the galaxy morphology, a natural predic¬ 
tion is that the number of satellites surrounding the galaxies should 
be also related with the shape of their host galaxies. It is worth not¬ 
ing that the merging activity is likely linked to the amount of mass 
contained within the bulge of the galaxies (e.g.|Hgmguist_&^aj 3 iej 


199lUSteinmetz & Muller 1995: Fu et al.l2003l:IZavala et all2008t 


Kroupa et alj 20ldl : iKroupal 20121) . The prominence of the bulge 
is one of the key ingredients in the galaxy morphological criteria, 
and consequently, a detailed segregation among the morphological 
types (beyond a disc/elliptical separation) could be connected with 
the number of satellites around the host galaxies. 

To classify our galaxies we have followed the traditional Hub¬ 
ble classification scheme. For all our galaxies, we look in detail the 
colour stamps provided for each of them by the finding chart tool 
of SDSS. To disentangle among ellipticals and SOs, we search for 
any evidence of a less steeply declining brightness in the outer re¬ 
gion of the galaxies beyond the central brightness condensation. In 
cases where the inclination of the disc component of the SO galaxy 
was clearly showing a flat outer structure or when some dust fea¬ 
tures were obvious, the distinction between the two galaxy types 
was relatively simple. Among the disc galaxy population, the seg¬ 
regation between Sa and Sb/c was done according to the relevance 
of the bulge in producing the overall light distribution as well as ex¬ 
ploring the properties of the spiral arm structure, i.e. the tightness 
with which the spiral arms are wound and the number of substruc¬ 
ture visible in those features. 


2.3 The sample of satellite galaxies 

To select our sample of potential satellites, we also use the 
‘specphoto’ catalogue of SDSS DR10. Our satellites are those 
galaxies in the catalogue which fulfil the proximity criteria (i.e. 
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projected distance to the hosts) and stellar mass criteria that we 
will explain in the next section (Section [3j. 

On building our sample of satellites, we have to prevent the in¬ 
clusion of objects with deficient measurement of its colour because 
this would lead to a wrong estimate of its stellar mass. To conduct 
this task, we need to account for the photometric errors both in g 
and r bands to assure the colour is measured with enough confi¬ 
dence. For this reason, in addition to the magnitude limit in the r 
band we have used above, we also demand that the photometric er¬ 
ror at estimating the number counts of each galaxy will be less than 
5cr the expected error at measuring their number counts. In other 
words, acceptable photometric error for each object for us are those 
whose error(counts) is <5x (counts + cr; ky ), with cr* y the uncer¬ 
tainty (in counts) at measuring the sky value in each band. We find 
as typical values for the sky in the SDSS images 24.88 counts (g 
band) and 23.96 counts (r band). We have used the following set of 
equations to transform our magnitudes and error(mag) provided by 
the catalogues into counts and error(counts): 



mag = 


2 5 log | counts | Q0.4(aa+Hxairmass) | 


exptime 


error(mag) = 


2.5 error(counts) 
In 10 counts 


(3) 

(4) 


with exptime=53.907 s and aa, kk and airmass provided for 
each object, being aa and kk the values of the zero-point and the 
extinction coefficient respectively. Those galaxies in our catalogue 
which show a photometric error larger than those values (in any 
of the two bands) are discarded from the analysis since these ones 
could be linked to artefacts in the image, proximity to bright nearby 
companions, etc. We have estimated the number of galaxies re¬ 
jected because of large photometric errors, finding that less than 
0.5 per cent of the objects are discarded, a reasonable result since 
these objects are relatively bright. 


Figure 3. Left-hand panel: the points in grey correspond to stellar mass dis¬ 
tribution versus redshift for all the galaxies in the ‘specphoto’ SDSS DR10 
catalogue. The pink points represent our sample of massive host galaxies. 
The black points indicate those objects in the right-hand panel within the 
redshift (0.023<z<0.025) used to estimate the stellar mass completeness. 
Right-hand panel: stellar mass distribution for the spectroscopic catalogue 
in the redshift interval 0.023<z<0.025. The red dashed line is the com¬ 
pleteness limit ~8.9xlO s M© whereas the blue line represents the minimum 
stellar mass for the sample of massive galaxies (10 n MO). A conservative 
estimation for the stellar mass completeness of the spectroscopic sample is 
showed by the orange dashed line: ~2.6xlO 9 M0 (see text for details). 


On what follows, we will consider this value as the most conserva¬ 
tive mass completeness limit of our satellite galaxies. 


2.3.1 The completeness of the satellite sample 

Once the stellar masses of the galaxies of our catalogue are deter¬ 
mined, we can estimate down to which stellar mass the catalogue 
is complete. The stellar mass limit for completeness is a function 
of the redshift (see Fig.[3]l. To explore the degree of completeness 
of satellites down to M* ~10 9 Mq, we have used our most distant 
redshift interval 0.023<z<0.025. The peak on the mass distribution 
of the galaxies in the catalogue at z-0.024 is log(M*/MO)~ 8.95. 
If we now take into account that the minimum stellar mass that we 
have fixed for our hosts (log(M*/Af©)~l 1.0), we should be able to 
study with completeness satellites whose Msac/Mnost <: 0.01. 

However, it is worth noting that there is a potential bias to 
miss the oldest satellites at a fixed stellar mass. This is because the 
catalogue is complete in redshift down to a given apparent r-band 
magnitude (r —17.7). That value translates into the following ab¬ 
solute magnitude for the sample at z=0.025: M r =-17.65 mag. To 
transform this absolute magnitude into a stellar mass limit we need 
to have an estimation of the stellar M/L ratio of our satellites. To 
do this, we assume a conservative age for the less massive galax¬ 
ies of 10 Gyr. Using the MIUSCAT spect ral en ergy distributions 
SEDs developed by Vazd ekis et alj J 201 2| ) and Ricciardell i et all 
d2012h . a solar metallicity and a Kroupa IMF, the (M/L) r ratio for 
these objects is around 3. This translates into the following stellar 
mass: log(M*/MO)~9.41 (i.e. a mass ratio satellite-host of 1:40). 


3 SATELLITE SELECTION CRITERIA 

Our criteria to search for potential satellites around the host sample 
are based on the next three steps. 

(i) We detect all the galaxies in the "specphoto’ SDSS DR10 
catalogue which are within a projected radial distance to our cen¬ 
tral galaxies of R =300 kpc. We only consider those host galaxies 
when the area enclosed by the search radius of satellites is fully 
contained within the catalogue borders. As we mentioned above, 
27 hosts were discarded due to their proximity to the catalogue 
edge. Our adopted search radius of 300 kpc is a compromise be¬ 
tween having a large area for finding a significant number of satel¬ 
lite candidates gravitationally bound to our central massive galaxies 
but not as large as to be severely contaminated by background and 
foreground objects (see Section I3.lt . 

(ii) The absolute difference between satellite redshifts and 
the redshift of the central galaxies must be lower than 1000 
km s -1 . This value has been used before in the literature to se- 
lect gravitationally bo u nd satellites of massive galaxies ( s ee e.g . 
IWang & Whitel 120121 : iRuiz. Trujillo & Marmol-Oueraltdl 12014]) . 
The velocity distribution of the galaxies around the massive host 
galaxies selected that way is close to a Gaussian shape with a dis¬ 
persion of 300 km s -1 . Consequently, our criteria enclose the vast 
majority of satellites around the massive galaxies. 

(iii) The mass ratio between our host massive galaxy and the 
satellite should be above 1:100. 

Those objects which fulfil the above criteria are counted as 
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potential satellites of their hosts. The number of satellites observed 
around a host sample is defined as N 0 b s ■ Before showing our results, 
we will address the potential biases that can affect our counting of 
satellites around the massive hosts. 


3.1 Background correction 

Despite we have used spectroscopic redshift information to select 
our potential satellite galaxies, there is still a fraction of objects that 
satisfy all the above criteria but are not gravitationally bound to our 
massive galaxies. These objects are counted as satellites because 
the uncertainties on their redshift estimates include them within our 
searching redshift range. These foreground and background objects 
(hereafter we will use the term background to refer to both of them) 
constitute an important source of uncertainty in this kind of stud¬ 
ies. Consequently, it is key to estimate accurately the background 
contamination in order to statistically subtract its contribution from 
the number of galaxies hosting satellites. 

To estimate the typical number of background objects that 
contaminates our satellite samples, we have developed a set of 
simulations. The procedure consists on placing a number of mock 
galaxies (equal to the number of our host galaxies) randomly 
through the volume of the catalogue conserving their original val¬ 
ues in the stellar mass and redshift of the sample of massive galax¬ 
ies. Once we have placed our mock galaxies through the cata¬ 
logue, we count which number of these galaxies have fake satellites 
around them taking into account the criteria of stellar mass, redshift 
and distance explained in the above section. This procedure is re¬ 
peated 2000 times to have a robust estimate of the number of 'satel¬ 
lites’ around the mock host galaxies. We define this average number 
as Ns- Then, being Nobs the number of observed satellite galaxies 
around either of our host massive galaxies, we correct statistically 
its excess subtracting the number of satellites representative of the 
background (Ns), such as it is shown in the equation below (equa¬ 
tion [5j. By construction, N s is independent on the morphology of 
the host. 

N Sat,S = (Vobs - Ns (5) 

To take into account that the environmental density and the 
morphology are linked for massive nearby galaxies, we also com¬ 
pute a clustering correction as explained in the following section. 


3.2 Clustering correction 


At low redshift, the over density regions are specially populated by 
massive galaxies. It is worth therefore exploring whether our back¬ 
ground correction is representative of the contamination of sources 
surrounding our host galaxies or whether it is necessary to com¬ 
pute the excess of probability of finding ’satellites’ in these envi¬ 
ronments. We term this as clustering. 

Being the clustering an effect associated with the region sur¬ 
rounding the hosts, ideally one would like to measure its influ¬ 
ence as closer as possible to the host. In practice, this is done 
by measuring the amount of s atellite candidates in different an- 
nuli beyond our search radius |chen_e£al. 2006; Liu et al^ 20~i~li; 
Marmol-Oueralto et alj 120121 : iRuiz. Trujillo & Marmol-Oueraltdl 
2014). We denote N c as the number of ’satellite’ galaxies placed in 
these annuli which fulfil our selection criteria for each morpholog¬ 
ical host type. N c is a measurement of the background contamina¬ 
tion plus the excess over this background caused by the clustering. 
This method has the disadvantage, compared to the simulations that 


we have conducted above, that is statistically more uncertain since 
Nc can be measured only around our massive galaxies and their 
number is relatively small. 

To quantify the effect of the clustering, we count the number 
Nc of satellites observed in annuli between 500 and 600 kpc which 
fulfil our selection criteria. The size of each annuli has the same 
area that our main exploration area around the hosts [i.e. n (300 
kpc) 2 ]. As it was done in the above section, we subtract Nc to the 
number of observed satellites N 0 b s to correct for the statistic excess 
given by the clustering in the sample of observed satellites. In other 
words: 

A Sat , c = N obs - Nc (6) 

The radial range 500-600 kpc for determining the clus¬ 
tering is a compromise among having a local measurement of 
the environment around our massive host galaxies but being 
far away enough such as the probability of finding a gravita¬ 
tionally bounded satellite to our targeted galaxy will be low. 
The projected radial distance of 500 k pc is cho sen follow- 


ing many works in the literature (e.g. 

Sales & Lambas 

2004 

2005: Chen et all 20061: 

Bailin et al. 2008 

; Wang & White 

2012 


iRuiz. Trujillo & Marmol-Oueraltdl 120141) which have used only 
galaxies with radial distances lower than 500 kpc to define their 
sample of truly (i.e. bounded) satellite galaxies. 


4 RESULTS 

4.1 Cumulative number of satellites per galaxy host 

The abundance of satellites is quantified using the number of satel¬ 
lites per number of massive galaxies Ns at /NHost down to a mass ratio 
satellite-host of 1:100. The results are shown in the left-hand panel 
of the Fig. [4] and Table Q] In addition to explore our satellites in a 
search radius of 300 kpc, we repeat the same exercise using a 100 
kpc radius to compare with previous results in the literature. 

The N Sa t/N H ost values have been corrected from contaminants 
by subtracting to Mobs the quantities Ns and Nc found in the back¬ 
ground simulation and in our estimates of clustering, respectively. 
The uncertainties of Mobs and N c were estimated from bootstrap 
resamplings of host sample sets. As we can see in Table Q] the 
background contamination is very low since our work uses spec¬ 
troscopic redshifts. According to that table, the maximum contam¬ 
ination expected by the background is ~0.1 satellites per galaxy 
host. The ratios N s /N H ost and N c /N Host in Table |T] show that the 
background contamination is, as expected, quite independent on 
the morphological type. However, if we compare the typical num¬ 
ber of fake satellites due to clustering (Nc/Nn ost ), we see significant 
differences. By far, the host galaxies that are more likely affected 
by contaminants are the E-type. On one hand, we have to take into 
account that the massive Sb/c-types are not usually expected to be 
into the cluster’s core regions. Thus, if we compare both M s /M Host 
and Mc/Mhosi around the massive Sb/c-types, we find a factor of 6 
larger contamination due to clustering than due to background. In 
contrast, this ratio increases dramatically up to ~16 times for the 
elliptical hosts. SO and Sa-types show a ratio of ~7-8 and ~6-7, re¬ 
spectively, a factor similar to the one found around Sb/c-types. It 
is remarkable that the density of objects at assessing the clustering 
and the background contaminant around the Sb/c-types, although 
close, are not similar. All this indicates that these massive galaxies 
are immersed in an environment similar to the samples of SO and 
Sa and therefore, they are not completely isolated. 
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Figure 4. Left-hand panel: cumulative number of satellites per galaxy host for each morphological type versus the stellar mass ratio satellite-host down to 
1:100. The dark red lines correspond to centrals E type, red to SO type, dark blue to Sa type and soft blue to Sb/c type. The coloured areas represent the space 
between the abundance estimated after applying background and clustering corrections. The dashed lines correspond to the background correction and the 
continuous line to the clustering. The striped grey region indicates our more conservative measure of the completeness (~ 1:40) — see Section ixm Right-hand 
panel: differential number of satellites per galaxy host for each of the morphological type versus the stellar mass ratio satellite-host. The colour code associated 
with the lines is the same that in the left-hand panel. 


Focusing on the satellite abundances, we find that the number 
of satellites per galaxy host is 1.4-2.0 when we explore satellites 
down to a mass ratio 1:10 around massive elliptical galaxies. If we 
increase that range of mass ratio down to 1:100, we find 3.7-5.4. 
The other morphological types show less number of satellites than 
in the E-host case. The more extreme case is found at comparing 
with the massive Sb/c type, their number of satellites per galaxy 
host grows between 0.23-0.46 and 0.88-1.63 from 1:10 to 1:100. 
Consequently, the massive E hosts have ~5 times more satellites 
than the massive Sb/c down to 1:10. A difference which declines 
up to 3-4 times in the case 1:100. At comparing Nsat/Nnost down 
to 1:10 with the other samples of galaxies, we find that the SO and 
Sa types host ~2 and 3 times less satellites, respectively than the 
E types. At extending our search down to 1:100, that difference 
among the SO-Sa types and ellipticals barely change. When we re¬ 
strict our satellite search up to only 100 kpc, the difference in the 
number of satellites among the different morphological types re¬ 
mains similar. As expected, the total number of satellites decreases 
when comparing a search radius of 300 kpc to one of 100 kpc. The 
decreasing factor is 2.5, 2.9, 2.2 and 3.5 for the E, SO, Sa and Sb/c 
types, respectively. 


4.2 Differential number of satellites per galaxy host 

We have also explored different intervals of the mass ratio satellite- 
host. The intervals are defined as 1:1-1:2, 1:2-1:5, 1:5-1:10, 1:10- 
1:20, 1:20-1:50, 1:50-1:100 and the number of satellites per host is 
illustrated in the right-hand panel of the Fig. [4] and Table [2] This 
‘differential’ test allows us to compare how is distributed the popu¬ 
lation of satellites respect to their stellar mass. 

As the results presented above, the dependence of the abun¬ 
dance of satellites with the morphology of the host is also appre¬ 


ciable. In general, the massive ellipticals have more satellites, fol¬ 
lowed by SO, Sa and Sb/c types. 

As shown in Table [2] the number of fake satellites per host 
due to the background As/Anost remains independently on the mor¬ 
phological type whereas the contamination due to the clustering 
Ac/A^uost is prominently larger in the elliptical case and very simi¬ 
lar among the types SO, Sa and Sb/c. 

4.3 The amount of mass surrounding the galaxy hosts 

We have studied the amount of stellar mass accumulated by the 
satellites around our samples of massive galaxies. This is estimated 
down to 1:100 and as a function of their projected distances to the 
central galaxy. This cumulative stellar mass of the satellites is mea¬ 
sured summing the stellar mass of all the satellites down to 1:100 
in each interval of the search radius up to 300 kpc and then, sub¬ 
tracting the amount of stellar mass in the fake satellites from the 
background simulations. To apply the clustering correction on the 
amount of stellar mass enclosed by the satellites, we estimate the 
amount of mass in ’satellites’ in the interval 500-600 kpc using an 
annuli with the same projected area than the ones used to study this 
quantity (i.e. from n 50 2 to n 300 2 ). The results are illustrated in 
Fig. [5] 

E-type massive galaxies are surrounded by a factor of 2-5 
more stellar mass respect to the rest of morphological types. If we 
repeat this exercise using only satellites up to 100 kpc, the stellar 
mass accumulated for lenticulars and spirals is similar, however, it 
is still a factor of 2-6 lower than the mass collected by ellipticals. 
Among the early types, the ratio of mass surrounding ellipticals to 
the mass surrounding lenticulars is typically placed around a factor 
of 2 larger. The spiral types show a more gradual accumulation of 
stellar mass with radius. It is worth noting that the amount of stel¬ 
lar mass enclosed by the satellites of massive ellipticals is almost 
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Table 1. Cumulative number of satellites per host galaxy using the background (Wsacs/Mrost) and clustering corrections Wsat.c/Mfost) within a radial distance 
of 300 kpc and down to a mass ratio satellite-host 1:100. Afobs/Mtost is the observed number of satellites per host. /Vs/A'iinsi is the number of satellites per host 
in our background simulation and Nc/N^ost the number of satellites per host within 500-600 kpc in our clustering analysis. The background and clustering 
corrections are estimated for each morphological type as it is shown in the table. 


37Sat/TfHost 

Afobs/Affost 

Afe/A/kost 

Nc/N Host 

A^sat,s /Afrost 

A'sat.c/Affost 

E 

0.50-1.0 

0.40 ± 0.03 

0.006 ± 0.0001 

0.080 ± 0.008 

0.36 ± 0.03 

0.29 ± 0.03 

0.20-1.0 

1.23 ±0.04 

0.021 ±0.0002 

0.325 ± 0.022 

1.19 ±0.04 

0.88 ±0.02 

0.10-1.0 

1.94 ±0.04 

0.038 ± 0.0003 

0.647 ± 0.046 

2.02 ± 0.04 

1.40 ±0.05 

0.05-1.0 

3.14 ±0.08 

0.056 ± 0.0004 

0.985 ±0.061 

3.10 ±0.07 

2.17 ±0.05 

0.02-1.0 

4.44 ±0.11 

0.085 ± 0.0008 

1.373 ±0.083 

4.40 ±0.13 

3.10 ±0.08 

0.01-1.0 

5.51 ±0.16 

0.113 ±0.0010 

1.840 ±0.105 

5.41 ±0.16 

3.67 ±0.11 

SO 

0.50-1.0 

0.23 ± 0.02 

0.005 ± 0.0001 

0.032 ± 0.010 

0.24 ± 0.02 

0.21 ±0.02 

0.20-1.0 

0.64 ± 0.04 

0.020 ± 0.0003 

0.168 ±0.019 

0.63 ± 0.04 

0.47 ± 0.05 

0.10-1.0 

0.90 ± 0.02 

0.037 ± 0.0006 

0.329 ±0.019 

0.95 ± 0.02 

0.63 ± 0.06 

0.05-1.0 

1.40 ±0.07 

0.056 ± 0.0008 

0.346 ± 0.020 

1.47 ±0.07 

1.13 ±0.07 

0.02-1.0 

2.35 ±0.06 

0.084 ± 0.0012 

0.629 ± 0.040 

2.39 ± 0.06 

1.82 ±0.08 

0.01-1.0 

2.99 ± 0.09 

0.113 ±0.0016 

0.804 ± 0.043 

3.11 ±0.09 

2.32 ±0.12 

Sa 

0.50-1.0 

0.11 ±0.01 

0.006 ± 0.0001 

0.017 ± 0.003 

0.12 ±0.01 

0.11 ±0.01 

0.20-1.0 

0.41 ±0.03 

0.022 ± 0.0002 

0.164 ±0.014 

0.40 ± 0.03 

0.23 ± 0.02 

0.10-1.0 

0.67 ± 0.03 

0.039 ± 0.0003 

0.301 ± 0.033 

0.67 ± 0.03 

0.38 ± 0.02 

0.05-1.0 

0.97 ± 0.03 

0.057 ± 0.0004 

0.349 ± 0.037 

0.98 ± 0.03 

0.64 ± 0.03 

0.02-1.0 

1.45 ±0.02 

0.086 ± 0.0006 

0.492 ± 0.053 

1.52 ±0.02 

1.05 ±0.05 

0.01-1.0 

1.89 ±0.04 

0.115 ±0.0008 

0.662 ± 0.056 

1.88 ±0.04 

1.26 ±0.06 

Sb/c 

0.50-1.0 

0.13 ±0.03 

0.006 ± 0.0001 

0.041 ±0.008 

0.15 ±0.03 

0.10 ±0.03 

0.20-1.0 

0.26 ± 0.03 

0.022 ± 0.0003 

0.158 ±0.015 

0.29 ± 0.03 

0.12 ±0.03 

0.10-1.0 

0.42 ± 0.05 

0.039 ± 0.0005 

0.218 ±0.019 

0.46 ± 0.05 

0.23 ± 0.04 

0.05-1.0 

0.55 ±0.06 

0.057 ± 0.0007 

0.368 ±0.019 

0.63 ± 0.06 

0.25 ± 0.07 

0.02-1.0 

1.02 ±0.10 

0.086 ± 0.0009 

0.555 ± 0.035 

1.17 ± 0.10 

0.60 ±0.12 

0.01-1.0 

1.45 ±0.12 

0.114 ±0.0010 

0.655 ± 0.023 

1.63 ±0.12 

0.88 ±0.13 


as large as the mass of the host galaxies (i.e. ~10 u Mq). The er¬ 
ror bars in Fig. [5] are estimated using the contribution of Poisson 
errors based on the number of observed satellites and the standard 
deviation of the average number of fake satellites found in the sim¬ 
ulations. 


4.4 The merging channel of massive galaxies 

Using the previous distribution of satellites, we can speculate about 
the stellar mass which could be potentially transferred to the hosts 
due to satellite infall. Under the assumption that eventually, all the 
satellites surrounding our massive galaxies will infall into their 
massive hosts, we can estimate which satellites could, in the fu¬ 
ture, contribute most to a potential mass growth of the host galaxy. 
On what follows, we assume that all satellites, independently of 
their mass, will infall with the same speed on the central galaxy. 
Note, however, that this could be not necessary true, since it is the¬ 
oretically expected that most massive satellites will have shorter 
merging time-scale (e.g. Ijiang et al.l2014ll . 

To probe the most likely merger channel what we have done 
is the following. We have added all the stellar mass contained by 
the satellites within the intervals of stellar mass studied, then, we 


divide this quantity by the sum of the mass of all the host galaxies, 
such as it indicates the following equation: 

^Sat-bin 

2 ^Sat-bin,! 


2 ^HostJ 

j= 1 

The sum of all the mass in the host galaxies is a fixed quantity 
for our samples of hosts and their values are 2^° sl ^Hostj=(l 1-2, 
5.26, 10.4, 6.76 )x10 12 Mq for (E, SO, Sa, Sb/c) types, respectively. 
These numbers correspond to the following typical masses per 
galaxy host: 2?T M HostJ /AW=(1.3, 1.4, 1.3, 1.3) xl0 u M o . 

Our results are illustrated in Fig. [6] and Table [3] for each host 
sample, once corrected by the effects of background and clustering. 
Fig.[6]highlights the different way the stellar mass around the mas¬ 
sive galaxies builds up as a function of the morphological type. The 
merger channel is mainly produced by satellites down to 1:10 for 
all the host samples. The average total amount of stellar mass con¬ 
tained in the satellite population down to 1:10 compared to the total 
amount of stellar mass in the hosts to the different morphological 
types studied is 45.4±6.6, 24.7±7.5, 14.3±3.8, 10.8±3.7 per cent 
for E, SO, Sa, Sb/c types, respectively. Down to 1:100 these values 
increase up to 54.5± 15.1, 30.5±8.7, 16.9±4.3, 12.6±4.5 per cent of 
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Table 2. Differential number of satellites per host galaxy using the background (iVsat.s INnost) and clustering corrections (A^sat.c/Mtost) within a radial distance 
of 300 kpc and down to a mass ratio satellite-host 1:100. Aiobs/^Host is the observed number of satellites per host. Ns/Nnost is the number of satellites per host 
in our background simulation and IVc/Ahosi the number of satellites per host within 500-600 kpc in our clustering analysis. The background and clustering 
corrections are derived independently for each morphological type as it is shown in the table. 


MSat/Tfilost 

Aobs/AHost 

As/A H „ st 

Ac/IVHost 

Asat,s/AHost 

Asat,c/AHost 

E 

0.50-1.00 

0.40 ± 0.03 

0.006 ± 0.0001 

0.080 ± 0.008 

0.36 ± 0.03 

0.29 ± 0.03 

0.20-0.50 

0.83 ± 0.03 

0.015 ± 0.0002 

0.251 ±0.022 

0.83 ± 0.03 

0.59 ±0.03 

0.10-0.20 

0.78 ± 0.03 

0.017 ± 0.0002 

0.322 ± 0.032 

0.83 ± 0.03 

0.52 ± 0.03 

0.05-0.10 

1.13 ±0.03 

0.018 ±0.0002 

0.338 ± 0.020 

1.07 ±0.03 

0.77 ± 0.02 

0.02-0.05 

1.36 ±0.05 

0.029 ± 0.0003 

0.389 ± 0.026 

1.30 ±0.05 

0.93 ± 0.04 

0.01-0.02 

0.94 ± 0.02 

0.028 ± 0.0003 

0.467 ± 0.027 

1.01 ±0.02 

0.58 ± 0.04 

SO 

0.50-1.00 

0.23 ± 0.02 

0.005 ± 0.0001 

0.032 ± 0.010 

0.24 ± 0.02 

0.21 ±0.02 

0.20-0.50 

0.41 ± 0.03 

0.014 ± 0.0002 

0.137 ±0.013 

0.39 ± 0.03 

0.26 ± 0.04 

0.10-0.20 

0.30 ± 0.03 

0.017 ± 0.0003 

0.174 ±0.018 

0.32 ± 0.03 

0.16 ±0.02 

0.05-0.10 

0.46 ± 0.04 

0.019 ± 0.0003 

0.005 ± 0.003 

0.53 ± 0.04 

0.50 ± 0.04 

0.02-0.05 

0.91 ±0.05 

0.028 ± 0.0004 

0.255 ±0.021 

0.92 ± 0.05 

0.68 ± 0.08 

0.01-0.02 

0.63 ± 0.04 

0.029 ± 0.0004 

0.141 ±0.016 

0.71 ±0.04 

0.50 ± 0.03 

Sa 

0.50-1.00 

0.11 ±0.01 

0.006 ± 0.0001 

0.017 ±0.003 

0.12 ±0.01 

0.11 ±0.01 

0.20-0.50 

0.27 ± 0.02 

0.016 ± 0.0002 

0.147 ±0.014 

0.27 ± 0.02 

0.12 ±0.01 

0.10-0.20 

0.29 ± 0.02 

0.017 ± 0.0002 

0.137 ±0.023 

0.27 ± 0.02 

0.15 ±0.03 

0.05-0.10 

0.30 ±0.02 

0.018 ±0.0001 

0.048 ± 0.009 

0.31 ±0.02 

0.26 ± 0.02 

0.02-0.05 

0.57 ± 0.03 

0.029 ± 0.0002 

0.143 ±0.017 

0.54 ± 0.03 

0.41 ±0.03 

0.01-0.02 

0.40 ± 0.02 

0.029 ± 0.0002 

0.156 ±0.014 

0.36 ± 0.02 

0.21 ±0.02 

Sb/c 

0.50-1.00 

0.13 ±0.03 

0.006 ± 0.0001 

0.041 ±0.008 

0.15 ±0.03 

0.10 ±0.03 

0.20-0.50 

0.13 ±0.01 

0.016 ± 0.0002 

0.125 ±0.014 

0.13 ±0.01 

0.02 ±0.01 

0.10-0.20 

0.14 ±0.02 

0.017 ± 0.0003 

0.046 ± 0.009 

0.17 ±0.02 

0.12 ±0.01 

0.05-0.10 

0.15 ±0.02 

0.019 ± 0.0002 

0.150 ±0.015 

0.17 ±0.02 

0.02 ± 0.04 

0.02-0.05 

0.51 ±0.06 

0.028 ± 0.0002 

0.186 ±0.020 

0.54 ± 0.06 

0.35 ± 0.05 

0.01-0.02 

0.43 ± 0.03 

0.028 ± 0.0002 

0.153 ±0.016 

0.46 ± 0.03 

0.29 ± 0.03 


the total amount of mass contained in their hosts. These numbers 
indicate that the contribution of the low mass population of satel¬ 
lites 1:10-1:100 to the host mass is not significant compared to the 
mass ratio 1:1-1:10. We have assumed Poisson errors to estimate 
our error bars in the Fig. [6] 

At limiting our exploration up to 100 kpc, the above re¬ 
sults decrease notably. The satellites more massive than 1:10 are 
again the dominant mass contributor. The percentage which re¬ 
flects that contribution with respect to the total mass of the hosts 
is roughly 14.6±3.2, 2.4±1.3, 7.9±1.0, 3.0±2.1 per cent for E, 
SO, Sa, Sb/c types, respectively, corresponding to the mean of 
background-clustering results once applied the corrections. Down 
to 1:100, these percentages increase up to 17.2±3.9, 4.5±2.1, 
8.7±2.8, 3.6±2.4. 

As we pointed out before, if the satellites eventually infall into 
the host galaxies, the merger channel will be largely dominated by 
satellites with a mass ratio below 1:10 for all morphological types. 
For satellites up to 300 kpc (100 kpc), this corresponds to 67.2, 
68.4, 88.1 and 85.7 (85, 57, 91 and 83) per cent of the total mass 
enclosed by the satellites within that radial distance. The contri¬ 
bution of the most massive satellites is particularly important for 
spiral types, and slightly weaker for early types. Down to 1:10, the 
distribution of stellar mass of the satellites reveals that the massive 


E galaxies could have their main contributor between the relative 
masses 1:5-1:1 whereas the main contributor of the SO, Sa and Sb/c 
types is the most massive satellites (1:2-1:1). It is also worth notic¬ 
ing that elliptical galaxies can have up to ~4 times more mass in 
satellites of 1:10-1:100 compared to Sb/c types. 

If the theoretical expectations are correct, and the merger time- 
scales are shorter for the most massive satellites, this mass growth 
due to the larger satellites will be even more important than the 
result shown in Fig. [6] 


5 DISCUSSION 

The results presented in this paper attempt to establish a robust 
z ~0 reference for the study of the evolution of the abundance of 
satellites around massive galaxies with cosmic time and their de¬ 
pendence with host morphology. In addition, it exte nds the r esults 
obtained in our previous work, iRuiz, Trujillo & Marmol-Oueraltdl 
{mi. in which we studied the abundance of satellites around E- 
type hosts using the SDSS Data Release 7 (DR7). In this paper, we 
have completed that results, using a deeper spectroscopic sample 
and studying the satellite abundances also according to their host 
morphology. 
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Figure 6. The merging channel of the massive galaxies as a function of their morphology. The panels show the contribution (in percentage) of the satellite 
mass enclosed in each mass bin to the total mass confined by their hosts for each morphological type. The red solid line represents this quantity after correcting 
for the background contaminant (T's), and the green solid line after correcting by clustering (Tc). The black, red and green dashed lines show 'I'ohs. T'siin 
and 'T ( - u found in the observations, background simulations and clustering, respectively. The numbers over each bin correspond to the number of observed 
satellites within each mass interval. The vertical dashed area illustrates the region where the satellite incompleteness could play a role. 


We have found that the abundance of satellites turns out to be 
significantly dependent on the morphology of the hosts for galax¬ 
ies with similar stellar masses. The abundance of satellites is much 
higher around elliptical galaxies. The fact that Sb/c galaxies have 
fewer satellites around them probably helps to explain how they 
have maintained their disc-like structure during their lifetime. On 
the contrary, the large number of satellites around massive ellip¬ 
ticals could help to clarify the characteristic large envelopes sur¬ 
rounding these objects that are thought to be connected with inten¬ 
sive accretion. 


Several works have recently estimated the number of satel¬ 
lites in different redshift intervals around massive galaxies both 


201(1 Guo et al] 

201 lbj, 

2013: Lares et al. 

2011; Liu et al. 

2011; 

Nierenberg et al. 

20 1J; 

Ouilis & Trujillo 

20121: TaletaQ 

2012; 


Wang & White! 20121 : ?). The observational studies agree within 


the measurement uncertainties. However, the satellites in the 
simulations are, in general, more abundant. Some of these previous 
studies have shown that the number of satellites around massive 
galaxies of a given mass ratio has remained constant since at least 


z ~2 

( Jackson et al. 20101; Nierenberg et al. 2011, 2013; Man et al. 

201 i 

; Marmol-Oueralto et al. 


2012L 20131; Newman et al. 

2012 

; Ouilis & Trujilloi 2012 


Huertas-Comnanv et al. 20131; 


Ferreras et al 1120141 ). In this paper, less affected by incompleteness 


at low stellar mass than previous works, we readdress this question 
and explore how the theoretical expectations compare to the 
observational data in the local Universe. We do this for our four 
different morphological types. 


5.1 Abundance of satellites: comparison with previous works 

5.1.1 Observations 

In order to check the agreement of our results with the litera¬ 
ture, we do a direct comparison with the abundance of satellites 
found by IWang & Whitel J20 1 2l ). In their work, the authors stud¬ 
ied the number of satellites around a large sample of isolated mas¬ 
sive galaxies to probe the ACDM scenario. In their paper, they 
segregated the massive galaxies not using visual morphology but 
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Figure 5. Cumulative stellar mass enclosed by satellites down to 1 : 100 as 
a function of the projected radial distance up to 300 kpc. The background 
contaminant correction, as well as the clustering correction, is applied for 
all the morphological types. The dark red lines correspond to centrals E, red 
SO, dark blue Sa and soft blue Sb/c type. 


colours, and estimated the abundance of satellites according to the 
stellar mass of the satellites instead of the mass_ratio satellite-host. 

To make a direct comparison with lWang & White] l2012h . we 
also estimate the abundance of satellites according _to the stellar 
mass instead the mass ratio satellite-host such as [Wang & White! 
HH did. Given that the average stellar mass of our host sam¬ 
ple is log(M*/Mg)~ll.l. we have compared our numbers with 
the average value they get combining their results from the host 
mass bins 10.8<log(M*/Mg)<ll.l and 1 l.l<log(M*/Mg)<l 1.4 
(their green and blue lines, respectively, in their figs 7 and 8). 
We confront our results for E and SO types with the ones they 
found for their red hosts. Additionally, we explore the differ¬ 
ence between the abundance of satellites in their blue hosts 
with our sample of disc galaxies (Sa and Sb/c types). Thus, 
we obtain, for satellites with log(M*/Mg)=10, Ns a INn osx = 0.31- 
0.44 (ours) and 0.14-0.19 (Wang & White) and, for satellites 
with log(M*/Mg)=9.0, N s JN Host = 0.50-0.65 (ours) and 0.44-0.55 
(Wang & White). When we compare our disc types with their 
blue massive galaxies we find, for satellites with log(M*/Mg)=10, 
(Vsat/Al Host = 0.05-0.1 (ours) and 0.07-0.12 (Wang & White) and, 
for satellites with log(M*/Mg)=9.0, IVsat/Mtost =0.22-0.29 (ours) 
and 0.08-0.15 (Wang & White). Our results are significantly 
larger at comparing their red primaries with our early-type galax¬ 
ies for log(M*/Mg)=10 and closer for less massive satellites 
(log(M*/Mg)=9). This trend is not reproduced in the comparison 
of blue primaries versus late-type galaxies where we find a closer 
agreement. __ 

The difference with lWang & Whitt] ( l2012h is even larger if we 
consider a more extreme comparison: our massive elliptical galax¬ 
ies versus their red primaries. In this case, we obtain for satel¬ 
lites with log(M*/Mg)=10, N s JN Kost = 0.41-0.61 (ours) and 0.14- 
0.19 (Wang & White) and, for satellites with log(M*/Mg)=9.0, 
(Vsat/Mfo S t= 0.53-0.78 (ours) and 0.44-0.55 (Wang & White). This 
highlights the importance of morphology when we study the num¬ 
ber of satellites. Several reasons can contribute to the diffe rences 
found at comparing the red massive galaxies of IWang & Whitel 


Table 3. The merging channel of local massive galaxies for each morpho¬ 
logical type. The table shows the contribution (in per cent) of the stellar 
mass enclosed in each satellite mass bin to the total mass confined by their 
hosts. We show the observed fraction (in per cent) Tobs as well as their 
values once corrected from background Tsat.s and clustering '1'sai.c- 


Msat/Tfilost 

^Obs 
(per cent) 

^Sat.S 
(per cent) 

^Sat,C 

(per cent) 

E 

0.50-1.00 

23.54 ±4.30 

23.19 ±4.23 

13.50 ±2.47 

0.20-0.50 

26.18 ±3.15 

25.73 ±3.10 

10.75 ± 1.29 

0.10-0.20 

12.26 ± 1.48 

12.03 ± 1.45 

5.04 ± 0.61 

0.05-0.10 

7.82 ±0.83 

7.69 ±0.82 

3.38 ±0.36 

0.02-0.05 

4.21 ± 0.40 

4.12 ±0.40 

1.56 ±0.15 

0.01-0.02 

1.56 ±0.17 

1.52 ± 0.17 

0.50 ± 0.06 

Total 

75.58 ± 10.33 

74.27 ± 10.16 

34.73 ±4.93 

SO 

0.50-1.00 

15.42 ±5.14 

15.09 ± 5.03 

14.23 ±4.74 

0.20-0.50 

9.95 ± 2.57 

9.51 ±2.46 

5.12 ± 1.32 

0.10-0.20 

4.44 ± 1.28 

4.21 ± 1.21 

1.23 ±0.36 

0.05-0.10 

3.73 ± 0.83 

3.59 ±0.80 

2.79 ± 0.62 

0.02-0.05 

2.91 ± 0.49 

2.82 ± 0.48 

1.23 ±0.21 

0.01-0.02 

0.96 ±0.18 

0.92 ±0.18 

0.28 ± 0.05 

Total 

37.40 ± 10.50 

36.14 ± 10.16 

24.89 ±7.31 

Sa 

0.50-1.00 

8.76 ±2.77 

8.35 ± 2.64 

7.57 ± 2.39 

0.20-0.50 

8.26 ± 1.76 

7.79 ±1.66 

0.09 ± 0.02 

0.10-0.20 

3.75 ±0.80 

3.51 ±0.75 

1.30 ±0.28 

0.05-0.10 

1.96 ±0.39 

1.83 ±0.37 

0.15 ±0.03 

0.02-0.05 

1.88 ±0.28 

1.79 ±0.27 

1.05 ±0.16 

0.01-0.02 

0.48 ± 0.09 

0.44 ± 0.08 

0.00 ± 0.01 

Total 

25.09 ±6.10 

23.70 ±5.77 

10.17 ±2.88 

Sb/c 

0.50-1.00 

10.43 ± 3.69 

10.04 ±3.55 

5.81 ±2.05 

0.20-0.50 

3.69 ± 1.40 

3.20 ±1.21 

0.00 ± 0.39 

0.10-0.20 

2.38 ± 0.79 

2.15 ±0.72 

0.36 ±0.12 

0.05-0.10 

1.14 ±0.38 

1.00 ±0.33 

0.00 ±0.18 

0.02-0.05 

1.65 ±0.31 

1.56 ±0.30 

0.56 ±0.11 

0.01-0.02 

0.60 ±0.12 

0.56 ±0.11 

0.06 ± 0.01 

Total 

19.89 ±6.69 

18.51 ±6.22 

6.78 ± 2.86 


d2012i) with our early-type galaxies. The most important is prob¬ 
ably the environment where the samples are immersed since their 
study focus only on isolated galaxies whereas we use all the galax¬ 
ies. However, the abundance of satellites is some closer at extend¬ 
ing the comparison to the poorer mass satellites. This is likely due 
to our incompleteness around 10 9 Mg. In contrast, this effect is not 
so evident at comparing blue-massive versus our late-type galaxies 
since these galaxies more likely live in less crowded environment. 

From the above results we highlight two different character¬ 
istics. First, the isolation criteria used in previous works probably 
biases the global distribution of satellites around the massive galax¬ 
ies. And second, the determination of abundance of satellites using 
morphological criteria enhances the differences in the number of 
satellites compared those based on segregating the samples simply 
by colour. This implies a stronger connection between the galaxy 
morphology and their number of satellites larger than a potential 
connection between the satellites and the stellar population of their 
hosts. 

We can also c ompare o ur re sults with the recent work by 
iRuiz. Trujillo & Marmol-Oueraltdl d2014l ) on the satellite popula¬ 
tions around massive elliptical galaxies in the local universe. In 
that work, we studied the abundance of satellites within a radial 
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distance of 100 kpc and selected the host samples within the red- 
shift range 0.02-0.065. We found, after applying the contamination 
corrections, that IVsat/Miost was between 0.24 and 0.28 for satel¬ 
lites up to 1:10 using a spectroscopic catalogue and N &ra /Nuost 
0.84 and 0.98 when we extended our search up to 1:100 using 
a catalogue of photometric redshifts (‘photoc’; SDSS DR7). In 
our new host sample of massive ellipticals, we obtain higher val¬ 
ues for both ratios (0.36-0.43 for 1:10 and 1.12-1.31 for 1:100) 
when we restrict our search radius down to R=100 kpc. The differ¬ 
ence between the average stellar mass o f both samples is 0.1 dex, 
being the average of Ruiz, Trujillo & MarmoTQueraltol ( 120141) ’s 
sample larger. This difference then could be attributed to techni¬ 
cal characteristics as the effect of the ’tiling’ of SDSS at differ¬ 
ent redshifts. In fact, in our previous work, we estimated that we 
could approximately lose up to 25 per cent of our potential satel¬ 
lites due to fibre collisions, a quantity which could partially ex¬ 
plain that difference in the case 1:10. Also, it could be playing 
a role the samples sele cted in b oth wo rks. In fac t, the selection 
of E galaxies jlone injRaiiz. Tru jillo & Marmol-Oueraltdl 120141) 
is based on iNair & Abraham! |2010|) whereas the present classi¬ 
fication is based on our own morphological analysis. Being the 
galaxies in the pre sent samp le sig nific antly closer than those in 
iRuiz. Trujillo & Marmol-Oueraltd72Q14h . our ability to distinguish 
among E and SO could be higher. In this sense, it worth stressing 
that the abundance of satellites found around SO massive galaxies 
for 1:10 and 1:100 is 0.12-0.15 and 0.68-0.72, respectively. If we 
combine the E and SO results, for 1:10, we get ~0.24-0.29 and 0.9- 
1.05 for 1:100. This is also in good agree ment with our previous 
results ( iRuiz, Trujillo & Marmol-Queraltcfeouil) . 

Leaving aside the number of satellites found per host, we can 
now draw our attention to the fraction of massive galaxies with 
satellites found in other works at z= 0, a param eter easily evalu¬ 
able from our study. In this context, iLiu et al] d201 lh . exploring 
Milky Way-like galaxies, found that only 12 per cent of those ob¬ 
jects have at least a satellite within R =100 kpc down to 1:10 mass 
ratio (private communication). We have evaluated the fraction of 
massive host having satellites down to 1:10 in our disc-like hosts. 
We find 10-12 per cent once corrected from clustering and back¬ 
ground contaminant. At segregating between the Sa and Sb/c sam¬ 
ples, that fraction is 0.15 and 0 .07, respectively. Our results are in 
good agreement with lLiu et alj ItoTII) . 

Also, we can make a comparison with iMarmol-Oueralto et al.1 
i2012l) . In this work, the authors conducted a similar analysis 
to what we have done here but for galaxies at higher redshifts 
(0.2<z<2). We compare our numbers with the galaxies they classi¬ 
fied as spheroids in their lower redshift range (0.2<z<0.75). They 
found that the fraction of massive galaxies with satellites around 
the sample and down to 1:10 is 23-28 per cent. Our results for 
spheroids, E and SO massive galaxies, are approximately 28 and 
34 per cent, respectively. As we commented before, this compari¬ 
son depends largely on the classification or grouping of the massive 
galaxies. A sample of spheroids with many SO objects reduces the 
satellite abundance whereas a higher purity in the E-type sample 
increases it. We can also compare with the galaxies they classified 
as disc-like in the same redshift range, they find that the fraction of 
massive galaxies with at least a satellite was 5-9 per cent. A closer 
but lower value to our disc-like 10-13 per cent. However, as we 
mentioned before, at splitting into the morphological types Sa and 
Sb/c (we obtained 15 and 7 per cent, respectively) the agreement 
could depend again on the number of members of each morpholog¬ 
ical type. 

Other works by IChenl d2008i ). I Guo et aO 1201 lbl ). 


iNie renberg et ah l2012h or Kawinwanichakij etalJ l2014ll also 
searched for satellites around massive galaxies segregating their 
samples by colours, early- and late-type morphologies or even 
quiescent and star-forming galaxies. Some of the results of these 
works can be compared with our results due to the similarity in 
the satellite search criteria. Nevertheless, these authors computed 
the potential satellites as a fun ction of the visual or r-band 
magnitude contr ast satellite-host JChenl 120081 : 1 Guo et ahl 1201 lbl : 
INierenberg et all2012l ) instead the mass ratio. Therefore, in order 
to compare with them, we also made our search for satellites using 
the r-band contrast magnitude between the satellite and its host. 


In lGuo et all J20 1 lbl) . the authors studied the luminosity func¬ 
tion of satellite galaxies around isolated bright hosts (M, ~ -21.25) 
using SDSS galaxy samples. They scaled their search radius to 300 
kpc and found that the average number of satellites per host was 
0.1 for red hosts and 0.06 for blue hosts when they studied satel¬ 
lites for Am, ~2.5 (see their fig.9). In addition, they split their sam¬ 
ple into early and late types, finding similar results. At compar¬ 
ing with their late or blue hosts, the results are far from our abun¬ 
dance of satellites using that r-band magnitude contrast, we obtain 
Nsat/NHost ~0.15-0.21. Also, when we compare with our results for 
early-type galaxies, we obtain Ns at/NHost ~0 .41-0.60, an amount 
significantly higher than those of Guo et all 1201 lbl ). A possible 
explanation to this difference could be linked to their stellar mass 
distribution. To build their fig.9 they use host galaxies within the 
magnitude range -20.75<M,, <-21.75 whereas our host sample is 
around -21.5<M,. <-22.5. 


Interestingly, IChenl 1200 8l) al so studied the sate llites around 
bright isolated galaxies, but unlike [Guo et ali 1201 lbl) . the authors 
chose their host sample in the nearby Universe (z< 0.045). That 
sample was somewhat brighter than IGuo et af] d201 lbll ’s one (- 
20<M,. <-23) and they used a mildly shallower sample of satellites 
(Am r < 2). Under these conditions, they found an abundance of 
satellites per host of 0.3-0.4 and ~ 1.0 for blue a nd re d hosts , respec¬ 
tively. This abundance of satellites found by lChenl 12008l ) is some 
lower than our results for Sa and Sb/c types combined 0.42-0.60 
for Am, <2 but close to our results around the Sb/c types (0.37- 
0.53) considering our clustering correction. Also, our abundance 
of satellites for (E, SO) hosts combined (1.31-1.79) is larger than 
which IcherJ l2008i) found. This likely due to the isolation criteria 
applied to select their host sample. 


As we have seen above, there are multiple factors which af¬ 
fect to comparison with other works. The usage of colours or mor¬ 
phology for selecting the hosts, the search criteria for satellites, the 
difference among the stellar mass distributions of confronted sam¬ 
ples or the density of galaxies where these samples are immersed. 
Nonetheless and in general, we get a reasonable agreement with 
the literature when we compare galaxies with disc like morpholo¬ 
gies, our Sa and Sb/c types. As we mentioned before, these mas¬ 
sive galaxies are thought to live in less crowded environments and 
therefore, their abundance of satellites should not be susceptible to 
change too much at applying them an isolation criterion unlike what 
happens with early-type hosts and mainly with the massive ellipti¬ 
cals. The majority of works showed above have used host samples 
using an isolation criterion. This probably biases the global distri¬ 
bution of satellites since it does not consider the excess of satellites 
linked to crowded regions in the real Universe. Precisely, in these 
comparisons, our abundance of satellites around early-type galax¬ 
ies is, typically, larger, especially considering our sample of mas¬ 
sive ellipticals. In fact, this excess is soft ened due to the mixture of 
E and SO galaxies iRuiz, Trujillo & Marmol-Queraltdl2014l) since 
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the SO massive galaxies are typically in similar environments than 
the Sa types. 

In that sense, a recent study bv lGuo et al.l d2015 ) investigated 
the dependence of the luminosity function of the host galaxies as 
a function of inhabiting a filament or not using the SDSS data set. 
They found that the filamentary environment can increase the abun¬ 
dance of the brightest satellites by a factor of ~2 compared with 
non-filament isolated galaxies. This can help to explain the dis- 
crepancy with other observational works which used isolation (e.g. 
IChenll2008l: iGuo et al.ll201 lbl : IWang & Whitell2012h since many of 
our massive elliptical galaxies are expected to be in dense regions. 

5.1.2 Theory 

lOuilis & Truiillcl d2012h estimated, using the Millennium Simu¬ 
lations, the expected fraction of massive galaxies with satellites 
with a mass ratio down to 1:10 and down to 1:100 within a 
sphere off? =100 kpc. They conducted their study exploring galax- 
ies from z=2 to now using three different semi-analytica l models 
dBower et alJ 20061 : IPeLucia et al.l2007l : lGuo et alj201 lah . At z=0, 
the theoretical expectations suggested that the fraction of massive 
galaxies with at least a satellite down to 1:10 ranges from 0.3 to 0.4 
and down to 1:100 from 0.6 to 0.7. We have studied that fraction 
in our work once applied the clustering corrections. E type (0.28, 
0.63), SO type (0.10,0.39), Sa type (044, 0.35), Sb/c type (0.06, 
0.20) for 1:10 and 1:100, respectively. For the full sample of mas¬ 
sive galaxies, we get 0.14 studying satellites up to 1:10 and 0.39 
down to 1:100. Summarizing, the theory makes a good prediction 
if we compare their numbers with massive ellipticals but it obtains 
an important over-prediction of satellites when we extend the host 
sample to all our host galaxies. 

It is worth stressing that this discrepancy among the theoreti¬ 
cal and the observational results can not be explained due to the dif¬ 
ferent volumes of exploration used in both works: a spatial sphere 
of 100 kpc in lOuilis & Truiillol d2012l) and a cylinder in redshift 
here. Because the way we have selected our galaxies in redshift, 
basically all the satellites in the line of sight of the host within a 
projected radial distance up to 100 kpc are taken. In that sense, 
at comparing with the Millennium Simulation our number of ob¬ 
served satellites should be an upper limit (as they are only restricted 
to 100 kpc in depth). As the number of theoretical satellites is larger 
than observed, we can confidently claim that there is a discrepancy 
with our observations. Fin ally, a potential loss of satellites in the 
work by lOuilis & Truiillol ( 120120 due to resolution effects will also 
increase the discrepancy between the simulations and the observa¬ 
tions. _ 

In this same theoretical context, we find the work of !Guo et al.1 
d2013i) who investigated the luminosity functions of galactic satel¬ 
lites around isolated bright hosts from model satellites placed into 
the Millennium and Millennium II dark matter simulations by the 
GALFORM semi-analytic galaxy formation model. These lumi¬ 
nosity funct ions allow us to compare with our results. In their work, 
IGuo et all d2013h used bright hosts within the range -21.5 < M r < 
-22.5 and the same search radius than us. They found that the num¬ 
ber of satellites per host was 0.29 and 0.07 (see the upper panel 
of their fig.8) segregating their sample in red and blue hosts for 
A m r ~2.5, respectively. Around this value, we find that the abun¬ 
dance of satellites is 0.15-0.21 for (Sa, Sb/c) galaxies. Our num¬ 
ber of satellit es per red host is also much larger than the found 
bv lGuo et alJ d2013l) . We obtain 0.41-0.60 for Am,. ~2.5. For less 
brighter satellites, A in r ~4, 1 Guo et aD d2013l) obtain for red and 
blue hosts 0.5 and 0.11, respectively, whereas we obtain the ranges 


0.40-0.72 and 0.18-0.29 for (E, SO) and (Sa, Sb/c), respectively. 
The abundance of satellites for our blue hosts re mains a factor of 2 
larger whereas the number of satellites found bv lGuo et alJ d2013l) 
increa se sig nificantly to reach a similar amount to ours. In their 
work, I Guo et aD d2013l ) also compared with observational results 
from SDSS DR8. Despite having found a good agreement between 
observations and models for red primaries, they found dramatic 
differences around blue primaries, placing the model a factor of 
2-3 fewer satellites than are present around comparable SDSS pri¬ 
maries. 

5.2 Mass and efficiency of the dark matter haloes 

Assuming that there is a link between the dark matter halo mass 
and the number of satellites a galaxy has |Wang & White 2 0121 : ?; 
iKawinwanic hakii et all |20 1 4| ). we can speculate how the relative 
abundance of satellites found in our different morphological sam¬ 
ples are related to their halo masses. For example, if we compare 
the combined abundance of satellites around E and SO types down 
to 1:100, with the number of satellites per host around Sb/c massive 
galaxies (grouping early versus late type), we obtain a factor of 3 
higher in the number of satellites around early-type galaxies than 
late types. This implies that on average, the dark matter halo mass 
of early massive galaxies could be three times larger than the halo 
masses associated with our Sb/c galaxies. We assume here that the 
number of sat ellites is approximatel y p roportional to the dar k mat¬ 
ter halo mass l lWang & White!20121 : ?). IWang & Whitel d2012l) also 
studied that difference using massive galaxies (log(M*/MQ)~l 1.2). 
In their work, they found that the red centrals had about a factor of 
2 more satellites th an blue centrals. _ 

Interestingly, iMandelbaum et all d20061) estimated the effi¬ 
ciency with which baryons in the halo of the galaxies have been 
converted into stars using a large sample of weak gravitational 
lenses (0.02<z<0.3), finding that the relative efficiency between 
late- and early-type galaxies with stellar masses (see their table 
3) between log(M*/Mo)~l 1.0 and 11.3 (the mean stellar mass of 
our samples is log(M*/Mo)~ll.ll) was 2.5-4.36. If we consider 
a similar stellar mass for our samples and using the Equation 7 of 
I Mandelbaum et all ( I2006I) . we can compare our relative conversion 
efficiency. That range of relative efficiency between late and early 
types is, considering our background and clustering corrections, 2.6 
and 2.9, a value in good agreement with their results. Following the 
above discussion, we can claim that the haloes of early-type mas¬ 
sive galaxies are typically 2-3 times less efficient to convert the 
baryons into stars than the haloes of late-type massive galaxies. 
A factor which rises, in our more extreme comparison, to 3-4 at 
comparing the sample of ellipticals and the Sb/c types. At extend¬ 
ing that relative efficiency in the terms explained before, to other 
morphological types, we find that our massive late spirals (Sa) are 
15-43 per cent less efficient than the Sb/c ones. Within early-type 
ones, the massive lenticulars are a 58-74 per cent more efficient 
producing stars than ellipticals ones. 

It is also interesting to compa re with the statistical approach 
carried out by iMoster et al.l d2010l) to determine the relationship 
between the stellar masses of galaxies and the masses of the dark 
matter haloes in which they reside. In their work, they estimated 
the average number of satellites as a function of the dark matter 
halo mass using a halo occupation model. The mean number of 
satellites is illustrated in their central panel of their fig. 10. Under 
the assumption that down to 300 kpc we are taking all the satellites 
around our massive galaxies and, considering that our completeness 
is around log(M/MO)~9.0 ( pointed line in their fig. 10), we could 
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do a direct comparison to establish the typical halo mass associ¬ 
ated with each one of our samples using the abundance of satel¬ 
lites obtained after the corrections of background and clustering. 
We obtain for E-type massive galaxies a halo mass M Ha io ~8.5- 
14.9x 10 I2 M o , for the SO types M Ha io ~5.4-7.2xl0 12 M o , for the 
Sa types M Ha io ~3.2-4.5x10 12 Mq and for Sb/c types a halo mass 
M Ha io ~2.2-3.5x10 12 Mq. 


From these data, we can conduct other interesting estimation 
using the average stellar mass jjf our samples of mass ive galax¬ 
ies and the equation (7) from iMandelbaum et all ( 20061 ). This al¬ 
lows us to estimate again the conversion efficiency of baryons 
rj obtained for the differen t morp hological types E, SO, Sa and 
Sb/c starting from the lMoster et alj J20 lfj )~s results. We find 0.06- 
0.10, 0.13-0.17, 0.19-0.27 and 0.24-0.39, respectively. Interest¬ 
ingly, at combining (E,S0) and (Sa,Sb/c), grouping in early and 
late-types, we find a relative conversion efficiency of 2.3-2.4 be¬ 
tween late and early-type galaxies. A factor close to the one 
already showe d before when we compa re that conversion effi¬ 
ciency with the Mandelbaum et al] d20061 )'s work. In this context, 
iMore et all j20 1 lh find that the difference between halo mass of 
red and blue centrals is ~0.4 dex as the stellar mass of the cen- 
tral is log(M < /Mfr ) )~ll.l. T his difference is also co mputed by 
Kawinwanichakij et~akl l2014h an d lPhillips et alj J20I 4|). They find 
a lower factor studying samples of quiescent and star-forming mas¬ 
sive galaxies with log(M*/MQ)> 10.78 (l<z<3) and log(M*/MQ)~ 
10.5 (locally). Their samples of quiescent centrals have a higher 
median halo mass by a factor of ~0.3 dex (factor 2). This lower 
factor compared to our results and the previous results found in the 
literature could be produced by the combination of two factors, a 
lower mean stellar mass of the galaxies explored, and a larger dif¬ 
ference among the halo masses of the samples studied in the nearby 
Universe respect to th e halo ma sses of the galaxies in the redshift 
range window of Kawinwa njchakii et al] a20 1 4I )’s sample. 

Finally, iDutton et al. l l201Cf ) estimated the star formation ef¬ 
ficiencies from satellite kinematics, weak gravitational lensing, 
and halo abundance matching at redshift z~0. They found that 
the formation efficiency of early-type galaxies reached a peak of 
~12 per cent at M* ~10 10 ' 5 It 2 Mq, decreasing to 2.8 per cent at 
M* ~10 1l4 1t 2 Mq. In contrast, this efficiency was between 26 and 
33 per cent for late-type galaxies whose stellar mean mass esti¬ 
mated was M* ~ 10 u '°tr 2 MQ. Both results are consistent with 
those ones showed before and therefore reinforce our hypothesis 
of proportionality between number of satellites and dark mass halo 
we assume as a starting point. 


5.3 The main contributor to the growth of massive galaxies 


There is growing consensus (see e.g. a discussion in lTruiillo et al .1 
that the size evolution can not be entirely explained by 
internal mechanisms l ike active galactic nuclei (AGN) activity 
dFan et al]|2008l . l20ld : iRasone-Figueroa & Granatdl20lih . How¬ 


ever, it is not clear what is the relevance of major versus minor 
merging in the growth of the galaxies. On one hand, major mergers 


(e.g. Ciotti & van Albada 200 1|; 

Nipoti. Londrillo. & Ciotti 20031: 

Bovlan-Kolchin, Ma, & Quataert 

2006; Naab et al. 2007) seem to 

be very 

scarce (at least since z~ 1: Bundv et al. 

2009; deRavel et al. 

2009; 

Wild et alj 2009; Lopez-Sanjuan et al. 

20101; Kaviraj et al. 


201 lh to play a major role in the growth of the galaxies. On the 
other hand, minor merging (favoured theor etically for its effi ciency 
on increasing the s ize of the galaxies: iKhochfar & Burkert 1 200a : 
iMaller et al.lI 2 OO 6 I : iHonkins et all l2009bl: iNaab et al] 120091) con- 


fron ts som e probl ems with the number of satellites found at z~l 
(e.g. lFerreras et al]|2014l) . 

In order to investigate the physical origi n of the observed 
strong increase in galaxy sizes since redshift z~2. IOseret ~aP) t2012h 
led a theoretical study in which they found that the evolution of 
massive early-type galaxies and their present-day properties are 
predominantly determined by frequent mergers of moderate mass 
(1:5) and not only by major mergers. iFerreras et iDioH probed 
the merging channel of massive galaxies (M* >10 n MO) over the 
z=0.3-1.3 redshift window and down to a mass ratio satellite-host 
1:100, segregating their sample into early and late-type massive 
galaxies (see their fig.9). They found that the main contributor to 
the growth of the host mass is those satellites whose mass ratio 
satellite-host is —1:3 for both samples. 

We find at z=0 a merger channel dominated by satellites with 
mass ratio satellite-host larger than 1:5 for E-type massive galax¬ 
ies, or dominated by satellites more massive than 1:2 if we consider 
our clustering correction as the most reliable. We find same result 
for the rest of morphological types, the most massive satellites are 
the main mass growth contributors. The growth of massive SO, Sa 
and Sb/c galaxies seems to be dominated by satellites more massive 
than 1:2. However, it is evident that massive ellipticals have a sig¬ 
nificant large number of satellites with masses ranging between 1:2 
and 1:5. Concretely, and if we focus on our results after clustering 
correction, the contribution of these satellites is similar to the ones 
within the range 1:1-1:2. This is not seen in the rest of morpholog¬ 
ical types. The SO types also show a non-negligible contribution to 
the host mass within the interval 1:2-1:5 but lower than a factor of 2 
compared to the contribution of satellites of similar mass around el¬ 
lipticals. In contrast, spiral types barely have satellites in this mass 
range. It is then interesting to study if this 1:2-1:5 merging channel 
remains growing when we assess the contribution of the satellites 
around a still more massive host sample of massive ellipticals and 
if it is reproduced at higher redshift. In this context, a better iden¬ 
tification of massive ellipticals and lenticulars may be key to check 
whether this highlighted channel 1:2-1:5 also exists. 

As the merger channel since z~ 1 is similar to the one found 
here locally, the observations suggest that the mass and size in¬ 
crease of the elliptical massive galaxies will be d omina ted by satel¬ 
lites with mass ratio within 1:1-1:5 (see alsolLopez-Saniuan et al.l 
l2012l : lRuiz. Trujillo & Marmol-Oueraltcil2014r) whereas other mor¬ 
phological types as SO, Sa and Sb/c seem to associate that increase 
to mergers with satellites with similar masses to the host. Note, 
however, that these statements assume that the merger time-scale 
are independent of the mass ratio between t he satellites a nd the 
host galaxies. More realistic scenarios (e.g. ljiang et ali20l4i) sug¬ 
gest that the merger time-scale rises as the mass ratio between both 
galaxies increases. Accordingly, the smaller satellites will take sig¬ 
nificantly more time to merge with the host galaxies than the more 
massive ones. 

Based on this, what we can claim with some confidence is that 
low-mass satellites with mass ratio below 1:10 would play a minor 
role in the mass increase of the host galaxies. They would be just 
very small in number to contribute to the mass growth, plus they 
will have very large time-scales to efficiently infall into the massive 
galaxies. However, the small satellites could be playing a major 
role in the construction of the stellar haloes of the galaxies (see e.g. 
ICooper et alj|2013h . If the theoretical expectation remains, and the 
merger time-scales are shorter for the most massive satellites, this 
mass growth due to the larger satellites will be even more important 
than the result showed in Fig. [6] 
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6 SUMMARY AND CONCLUSIONS 

In this paper we explore the abundance of satellites around 254 
massive (10 11 <M* <2xlO n MO) low-;; (z < 0.025) galaxies visu¬ 
ally classified as E, SO, Sa and Sb/c. Using the SDSS DR10 spectro¬ 
scopic catalogue, the proximity of our host galaxies guarantees that 
we can explore satellites with completeness down to M* ~ 10 9 M©. 
Our satellite galaxies have been identified within a projected radial 
distance of 300 kpc around the central galaxy. A careful statistical 
analysis of the background and clustering has been applied to de¬ 
contaminate the number of satellites from fake satellites. The abun¬ 
dance of satellites decline significantly from the E galaxies to SO, 
Sa and Sb/c types showing an important dependence with the mor¬ 
phology of the host independently of the mass ratio satellite-host. 
The average number of satellites down to a mass ratio 1:100 within 
300 kpc is 4.5 ±0.3 for E hosts, 2.6 ±0.2 for SO, 1.5 ±0.1 for Sa and 
1.2 ± 0.2 for Sb/c. These quantities decrease by a factor of 2.5-3.5 
down to 1:10. 

Under the assumption that there is a proportionality between 
the number of satellites found and the dark matter halo mass, we 
find that the haloes of massive ellipticals are less efficient than the 
haloes of Sb/c types on converting their baryons into stars by a 
factor of 3-4. We need a halo 3-4 times more massive to create an 
elliptical with the same stellar mass than a Sb/c spiral. This factor 
decreases to ~3 when we group our samples in early- and late-type 
massive galaxies. 

If the satellites would eventually infall into their host galaxies, 
the growth of massive galaxies will be dominated by satellites with 
a mass ratio down to 1:10. Those satellites are the main contributors 
to the stellar mass enclosed by the satellites and responsible of the 
67.2, 68.4, 88.1 and 85.7 per cent of the total mass in satellites for 
E, SO, Sa, and Sb/c types, respectively, down to 300 kpc. Massive 
ellipticals seem to be surrounded by a remarkably larger number 
of poor massive satellites whereas the rest of morphological types 
typically merge with more massive objects. Specifically, the main 
contributor to the growth of massive spirals seems to be the satel¬ 
lites more massive than 1:2 (to the SO, Sa and Sb/c types). To the E 
hosts, the merger c hannel peaks wit hin 1:1-1:5, in agreement with 
the 1:5 pointed by lOser et alJ d2012h . 

These results could be used in future works to test the ACDM 
predictions about the number of satellites surrounding the most 
massive galaxies in the present-day Universe according to their 
morphology. In addition, the results presented here show the most 
likely merging channel of present-day massive galaxies. Finally, 
our work highlights the importance of the environment where mas¬ 
sive galaxies are immersed and how that environment is strongly 
linked to the host morphology. 
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